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SECTION  I 


INTRODUCTION 

In  air  defense  and  air  traffic  control  systems,  data  from  the 
system  radars  are  stereographically  prolected  onto  a common  coordi- 
nate plane  for  presentation  to  system  operators.  The  stereographic 
projection  of  radar  data  involves  two  steps;  stereographic  projection 
vising  slant  range,  azimuth  and  height  information  to  obtain  polar 
coordinates  iu  a plane  of  projection  centered  at  the  radar  site,  and 
transformation  of  the  radar  coordinates  into  cartesian  coordinates 
on  a common  coordinate  plane. 

The  stereographic  projection  and  transformation  process  is 
mathematically  complex.  Because  of  the  computational  complexity, 
several  assumptions  and  approximations  have  been  made  to  expedite 
processing  time  without  unduly  sacrificing  accuracy.  The  present 
SAC.E/BOIC  equations  for  computing  radar  coordinates,  although 
satisfactory  for  their  intended  usage,  introduce  unacceptable  regis- 
tration errors  when  extended  to  large  regions  as  will  be  encountered 
in  the  Joint  Surveillance  System  (JSS). 

11»is  report  describes  analysis  that  was  performed  on  the  stereo- 
graphic  projection  process.  Equations  for  obtaining  stereographic 
ground  range  are  derived.  The  derivation  indicates  that  the  SAGE/ 
BUIC  ground  range  equation  lacks  a scale  factor  vital  to  proper 
registration  in  large  regions.  The  scale  factor  is  a function  of 
the  radius  of  the  earth  at  a radar  site  and  the  radius  of  the  con- 
formal sphere. 
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EFFECT  OF  THE  CONFORMAL  SPHERE  ON  STEREOGRAPHIC  PROJECTION 


STEREOGRAPHIC  PROJECTION 


Stereographic  projection  is  a method  for  mapping  points  in 
space  onto  a plane  tangent  to  a sphere.  This  sphere  is  termed  the 
conformal  sphere  and  its  radius  is  arbitrary.  Figure  I depicts  the 
projection  geometry  for  a cross  section  of  the  sphere.  The  cross 
section  is  obtained  by  passing  a plane  through  the  center  of  the 
sphere,  the  aircraft,  and  the  point  of  tangencv . The  intersection 
of  this  plane  and  the  sphere  is  a great  circle.  Mapping  of  point  A 
in  space  onto  the  tangent  plane  BB'  results  in  the  stereographic 
ground  range  DT.  DT  is  obtained  by  passing  a line  from  point  0, 
opposite  the  point  of  tangency  D,  through  point  P,  the  point  of 
projection,  and  intersecting  the  tangent  plane.  Point  P is  the 
intersection  of  the  line  containing  point  A and  the  center  of  the 
earth,  point  C,  and  the  great  circle.  If  angle  DCP  is  designated  as 
ij>,  then  angle  DOP  equals  iJ>/2  since  DC=C0=CP.  The  stereographic 
ground  range  DT  is  determined  as  follows . 

DT  - DO  tan  (1) 


If  DT  is  defined  as  R and  DC,  CO,  CP  are  defined  as  Ec>  then 
equation  (1)  takes  the  form: 


R 


2Ec  ta"  (I) 


s 2E 

c 


1 - cosijj 
1 + cosij; 


■1/2 

. 


(2) 


The  stereographic  ground  range  R is  therefore  directly  proportional 

to  the  radius  of  the  conformal  sphere  E . 

c 
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THE  RADAR  PLANE 


Radar  data  are  stereographically  projected  onto  a plane  centered 
at  the  reporting  radar  site  and  tangent  to  the  conformal  sphere.  The 
plane  will  be  termed  the  radar  plane.  The  coordinate  axes  of  the 
radar  plane  are  oriented  such  that  the  positive  Y axis  is  directed 
toward  true  north  and  the  positive  X axis  towards  east.  Figure  II 
shows  the  orientation  of  the  coordinate  axes  on  the  radar  plane.  An 
aircraft's  location  on  the  plane  is  expressed  in  polar  coordinates. 

The  range  R is  the  stereographic  ground  range,  and  the  azimuth  angle 
9 is  the  azimuth  of  the  radar  return.  Azimuth  angles  are  measured 
clockwise  from  the  positive  Y axis. 

DETERMINATION  OF  STEREOGRAPHIC  GROUND  RANGE 

Radar  slant  range,  height  of  the  aircraft  above  sea  level  and 
elevation  of  the  site  above  sea  level  are  used  to  determine  the 
stereographic  ground  range  of  an  aircraft  on  the  radar  plane.  The 
angle  <p  between  the  radar  site,  the  center  of  the  earth  and  the 
aircraft  are  used  in  equation  (2)  to  calculate  the  stereographic 
ground  range.  The  angle  i|i  can  be  calculated  if  the  earth  is  assumed 
to  be  spherical.  Figure  III  illustrates  the  radar  geometry  for  a 
cross  section  of  a spherical  earth  where  Eg  is  the  radius  of  the 
earth,  h is  the  site  elevation,  S is  the  measured  slant  range  and  H 
is  the  aircraft  height.  From  Figure  III  the  angle  may  be  calculated 
from  the  law  of  cosines. 

S2  » (E  + h)2  + 
s 

cos4>  ■ 


(E  + H)  - 2(E  + h)(E  + H)  cosiji 
s s s 

2 2 

1 4.  r_h)__-  s"  __ 

2(E  -f  h)  (E  + H) 

S 8 


l - 


2(E  + b)(K  >•  H) 

b 8 


(3) 
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where: 


2 2 2 
F - S - (H  - h) 


Substituting  the  results  of  equation  (3)  into  equation  (2), 
the  ground  range  R may  be  obtained  as  follows. 

r.2  ll/2 


R - 2E 


E F 
c 


4(E  + h)(E  + H)  - F 
s s 


1 4.  H 4-  h Hh 

E r2 

9 E 

s 


4E 


-1/2 


(4) 


Reference  1 page  7 presents  the  SAGE/BUIC  formulation  of  the  ground 
range  wherein  it  is  assumed  that  Ec  and  Eg  are  equal  and  thus  cancel.  It 
will  be  shown  that  for  the  large  regions  that  will  be  encountered  in  the 
JSS  system,  this  assumption  results  in  unacceptable  registration  errors. 

APPROXIMATIONS  TO  GROUND  RANGE 

The  accuracy  of  typical  common  digitizer  search  radar  out- 
puts is  0.25  nmi  in  range  and  0.18°  in  azimuth.  Registration 
errors  are  a combination  of  data  errors,  radar  site  location  errors, 
and  errors  due  to  approximations  in  the  stereographic  projection 
process.  Since  equation  (4)  is  computationally  complex,  an  approxi- 
mation which  does  not  unduly  sacrifice  accuracy  is  used  to  expedite 
processing  time.  A maximum  error,  induced  by  approximation, of  0.18 
nmi  provides  a reasonable  compromise  between  processing  requirements 
and  registration  accuracy.  Four  approximations  are  presented  in  this 
section;  the  series  approximation,  the  first  order  approximation, 
the  JSS  approximation  and  the  current  SAGE/BUIC  approximation. 
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Equation  (4)  may  be  written  in  the.  following  form. 


h P 


The  term  within  brackets  in  the  denominator  of  equation  (5)  may  then 
be  expressed  by  the  following  series  expansion. 


The  bracketed  term  in  the  demonlnator  of  equation  (6)  contains 

a first  order  terra  and  two  second  order  terms.  The  maximum  value  of 

— j--—  is  0.00248,  the  maximum  value  of  is  0.000000721,  and  the 

s 2E„ 

2 3 
F 

maximum  value  of  — x is  0.000424.  Since  the  second  order  terms  are 

8EZ 

s 

smaller  than  the  first  order  terra,  they  may  be  neglected  and  equation 
(6)  may  be  expressed  as  follows. 
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R = 


E F 
c 


(7) 


The  JSS  Approximation 

By  replacing  the  aircraft  height  terra  (H)  in  the  denoninator  of 
equation  (7)  by  a constant,  equation  (7)  nay  be  expressed  as  follows. 


R a 


E F 
c 


where: 


(8) 


Hg  Is  the  raaxiraura  expected  aircraft  altitude  equal  to  100,000  ft. 
Equation  (8)  may  be  expressed  In  the  following  form. 

R « CP 

where:  v 


C » 


F5^ 


This  la  a particularly  good  approximation  since  the  stereographic 
ground  range  la  obtained  from  a simple  scale  multiplication  of  the 
quantity  F.  This  (treat ly  decreases  the  time  reauired  to  process 
radar  returns.  The  constant  C la  adaptation  defined  on  a alte-by- 
alte  bails. 

The  Current  SAflB/BUIC  Approximation 

If  E is  assumed  to  be  the  same  as  E , equation  (8)  la  expressed 

C R 

as  follows. 


R * 


(■  •*£) 


(9) 
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EARTH  MODEL 


The  earth  is  not  a perfect  sphere.  Therefore*  for  precise 
calculations  of  stereographic  ground  range*  a model  for  the  geometric 
shape  of  the  earth  must  be  adopted.  An  appropriate  first  order 
representation  is  an  ellipsoid.  The  ellipsoid  is  generated  by 
revolving  an  ellipse  about  its  seaiainor  axis.  The  earth  model  can 
therefore  be  specified  by  its  sealaajor  axis  or  equatorial  redlus 
Eq  and  the  eccentricity  e.  A cross  section  of  the  adopted  earth 
model  is  shown  in  Figure  IV.  The  eccentricity  is  defined  as  follows. 

e2  - 2f  - f2  (10) 

where: 


E is  the  seminlnor  axis  or  polar  rad lux 
P 

The  International  Ellipsoid  of  1924  will  be  used  for  the  purpose  of 

2 

this  report;  thus*  equals  3444.054  mi  and  e equals  ,00622267. 

Happing  from  Ellipsoid  to  Sphere 

The  etercographic  projection  equations  have  been  derived  for  a 
sphere.  It  Is  therefore  necessary  to  transform  points  on  or  above 
the  ellipsoid  to  points  on  or  above  the  sphere.  This  transformation 
must  be  conformal,  i.e.,  angle  preserving*  if  the  final  stereographic 
projection  is  to  be  conformal.  Reference  2 page  86  derives  the  re- 
lationship between  the  ellipsoid  and  the  conformal  sphere.  The 
mapping  of  points  on  or  above  a location  on  the  ellipsoid  onto  the 
conformal  sphere  is  performed  as  follows. 


j 

i 
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Figure  IZ  A CROSS  SECTION  OF  THE  ADOPTED  EARTH  MODEL 


In  these  equations: 

L,  X are  the  geographical  latitude  and  longitude  of  the  point 
on  the  ellipsoid 

are  the  latitude  and  longitude  of  the  corresponding 
point  on  the  conformal  sphere.  The  latitude  0 is 
called  the  conformal  latitude, 
e is  the  eccentricity  of  the  earth. 

Since  the  mapping  process  involves  only  a transformation  of  latltudec, 
the  stereographic  projection  equations  are  valid  for  the  ellipsoid 
if  the  conformal  latitude  (4>)  is  used  in  place  of  the  geographic 
latitude  (L). 


Radius  of  the  Earth  to  a Radar  Site 

In  determination  of  the  stereographic  ground  range,  a spherical 
earth  was  used  to  compute  the  angle  that  subtends  the  radar  site 
and  the  target.  Since  the  earth  is  actually  modeled  bv  an  ellipsoid, 
the  use  of  a spherical  earth  in  calculating  stereographic  ground 
range  without  introducing  corrections  to  slant  range  and  height^  to 
allow  for  the  conformal  projection, will  introduce  a certain  auount 
of  error.  To  minimire  this  error,  the  radius  of  the  spherical  earth 
Eg  is  set  equal  to  the  distance  from  the  center  of  the  ellipsoid  to 
the  surface  of  the  ellipsoid  at  the  radar  site  as  shown  in  Figure 
IV.  The  distance  E is  calculated  as  follows. 

E - <X2  4 Y2)1/2  (12) 

s is 


where: 

X “ N cosL 

H 

T ■ (1  - e2)  N slnL 

3 
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N 


E 

(1  - e2sin\)l/2 
L is  Che  geographic  latitude. 

Table  I gives  the  radius  of  the  Earth  at  latitudes  between  0 and 
90  degrees. 


Table  X 

Earth  Radius  Versus  Latitude 


Geographic  Latitude 

Earth  Radius  (nni) 

0° 

3444.054 

10* 

3443.707 

20* 

3442.; 08 

30* 

3441.173 

40* 

3439,286 

50* 

3437,273 

60* 

3435.375 

70* 

3433.824 

80* 

3432,810 

90* 

3432.458 

Calculations  show  that 

the  resulting  error  is  greatest  for 

« target  north  or  south  of  the  radar  and  does  not  exceed  + 0.03  wai 
at  a range  of  230  nai. 

THE  COMMON  COORDINATE  PLANE 

tn  a large  air  surveillance  region  such  at*  will  be  encountered 
in  the  JSS,  several  radara  are  United  together  in  order  to  diaplay  a 
composite  air  surveillance  picture.  It  la  therefore  necessary  to 
transfers  coordinates  on  the  individual  raoar  planes  to  coordinates 
on  a cosaon  plana.  The  transform*  t ion  process  requites  that  the 
coordinates  of  the  rsdara  on  the  coamon  plane  and  that  angular  ro- 
tations between  the  radar  planes  and  Che  cotton  plane  be  known. 
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The  origin  of  the  common  plane  or  region  center  is  defined  as 
the  center  of  the  smallest  circle  that  will  circumscribe  all  radars 
tied  into  the  air  surveillance  region.  The  common  plane  is  tangent 
to  the  conformal  sphere.  The  coordinates  axes  are  oriented  such 
that  the  positive  Y axis  is  directed  towards  true  north  and  the 
positive  X axis  is  directed  towards  east.  Coordinates  on  the  coosaon 
plane  are  expressed  in  cartesian  coordinates. 


Radar  Coordinates  on  the  Common  Plane 


A radar  site  or  other  known  location  can  be  stereographicallv 
projected  onto  the  common  coordinate  plane.  Reference  3 page  53 
derives  the  equations  necessary  to  project  a point  on  the  conformal 
sphere  onto  the  common  coordinate  plane.  The  rectangular  coordinates 
of  a radar  site  or  other  known  location  Xr»  •¥  on  the  common  plane 
are  obtained  as  follows. 

V „ 2E  — _ . (13) 

r c 1 + sin$sin$Q  + cos$cos$0cosAA 


sinSeosA-  * cag^nin^cesA* 
y « 2E  — — — - -- — ■ — -it — — — — - 

r “cl  + sinisinj  + cos^cosO^cos&A 


(IV) 


where: 


G,  A 


V X0 


AA  - A0  - A 


■ A - A. 


are  the  conformal  latitude  and  longitude  of  the 
point  to  he  projected 

are  the  conformal  latitude  and  longitude  of  the 
region  center 

if  longitudes  are  measured  positive  west  of  the 
prime  meridian 

if  longitudes  are  measured  positive  east  of  the 
prime  Meridian 
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Angular  Rotation 


As  indicated  in  reference  3 page  6,  an  angular  rotation  of  the 
radar  plane  with  respect  to  the  common  coordinate  plane  is  necessary 
for  the  transformation  process.  The  effect  of  the  rotation  is  to 
make  the  axes  of  the  radar  plane  more  nearly  parallel  to  the  axes  of 
the  common  plane.  The  angle  of  rotation  is  shown  in  reference  3 to 


be: 


B * tan 


-1 


- (sln$  + ain$g)sinAA 
cosAcos$0  + (1  + sin$sin$Q)cosAA 


(15) 


Transformation  of  Radar  Coordinates  on  the  Common  Plane 

The  equations  for  transformation  of  radar  coordinates  to 
coordinates  in  the  common  coordinate  plane  are  derived  in  reference 
3 pages  3 - 15.  The  exact  transformation  equations  involve  an  in- 
finite series.  To  lessen  the  processing  requirement  without  unduly 
sacrificing  accuracy  a second  order  approximation  is  used.  Rectangular 
coordinates  X,  V are  obtained  as  follows. 


where: 


X • Xf  *f  X(Rsin(e  + 0)  + AR2sin(2(9  + 0)  - y]) 

V - Y ♦ K(fceos(0  + 0)  + AR2cO3[2(0  + 0)  - yj) 
r 


K • 1 + 


4E 


A * 


4B 


tt  - (X2  + Y2)i/2 


(16) 

(1?) 
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R is  the  stereographic  ground  range 

0 is  the  azimuth  angle  measured  clockwise  from  north  at  the 
radar  site 

Xr,  are  the  coordinates  of  the  reporting  radar  on  the 
common  plane 

Since  the  coordinates  of  an  aircraft  on  the  common  plane  are  a 
function  of  the  stereographic  ground  range  Rf  any  error  in  the 
stereographic  ground  range  will  appear  as  a misregistration  on  the 
common  coordinate  plane. 


THE  CONFORMAL  SPHERE 


The  common  coordinate  plane  and  all  radar  planes  are  tangent  to 
the  conformal  sphere.  The  radius  of  the  conformal  sphere  is  arbi- 
trary, but  is  chosen  to  minimize  the  scale  errors  that  will  be  en- 
countered in  the  air  surveillance  region.  Scale  errors  result  in 
t.  , . tng  from  the  ellipsoid  to  the  conformal  sphere  and  in  mapping 
from  the  co~formal  sphere  onto  the  tangent  plane. 


Scale  Factor  ~ Ellipsoid  to  Conformal  Sphere 


The  scale  factor  associated  with  mapping  from  the  ellipsoid  to 
the  conformal  sphere  is  the  ratio  of  a linear  element  on  the  conformal 
sphere  to  a corresponding  linear  element  on  the  ellipsoid.  From 
reference  2 page  86,  the  ^cale  factor  associated  with  a point  of 
projection  is: 


Eccos<f> 
^1  NcosL 


(18) 


where : 


<j>  is  the  conformal  latitude  of  the  point  to  be  prelected 
L is  the  geographic  latitude  of  the  point  to  be  projected 


The  scale  factor  is  a function  of  the  radius  of  the  conformal  sphere 
and  the  conformal  and  geographic  latitude  of  the  point  to  be  pro- 
jected. The  scale  factor  can  be.  expressed  as  follows. 


Ec  (1  - e28in\)^cos4 


q 

E 

- -£k» 
E 1 


cosL 


(19) 


Table  II  gives  values  of  k^  for  several  different  latitudes. 


Table  II 

Scale  Factor  ~ Ellipsoid  to  Sphere  versus  latitude 


Geographic  Latitude 

0° 

10° 

20° 

30° 

408 

50* 

60° 

70® 

80® 

90® 


Scale  Factor  k^ 

1.00000000 

1.00010071 

1.00039099 

1.00083667 

1.00138493 

1.00197024 

1.00252203 

1.00297314 

1.00326821 

1.00337838 


Scale  Factor  ~ Conformal  Sphere  to  Plane 

The  scale  factor  associated  with  mapping  from  the  conformal 
sphere  to  the  plane  of  projection  is  the  ratio  of  a linear  element 
on  the  plane  of  projection  to  a corresponding  linear  element  on  the 
conformal  sphere,  The  scale  factor  is  therefore  a function  of  the 
angular  separation  between  the  origin  of  the  plane  and  the  point  of 
projection.  The  distance  on  the  sphere  between  the  origin  of  the 
plane  and  the  point  of  projection  is  calculated  as  follows. 

D - Ec<|*  (20) 
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The  stereographic  ground  range  R is  given  by  equation  (1)  as 
follows . 


R - 2Ectan(|)  (21) 

The  scale  factor  associated  with  mapping  from  the  conformal  sphere 
to  the  plane  of  projection  is  calculated  as  follows. 


Table  III  presents  the  scale  factor  k£  as  a function  of  angular 
separation. 

Table  III 

Scale  Factor  — Conformal  Sphere  to  Plane  versus  Angular  Separation 


Angular  Separation 

Scale  factor  k^ 

0° 

1.00000000 

2° 

1.00030468 

4° 

1.00121946 

6° 

1.00274658 

8° 

1.00488976 

10° 

1.00765427 

: 12° 

1.01104690 

14° 

1.01507605 

! 16° 

1.01975173 

18° 

1.02508563 

20° 

1.03109120 

Aircraft  and  radar  locations  in  the  air  surveillance  region  are 
presented  on  the  common  plane.  The  scale  factor  kj,  associated  with 
the  common  plane  is  a function  of  the  angular  separation  between  the 
region  center  and  the  point  of  Interest.  From  reference  3 page  23, 
the  angular  separation  between  the  region  center  and  a point  of 

l 
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interest  can  be  expressed  in  terms  of  their  locations  on  the  conformal 
sphere  as  follows. 


cosifi  ■ simJisin^Q  4-  cos^cos^cosAX  (23) 

Combining  equations  (22)  and  (23)  the  scale  factor  is  expressed 
as  follows. 


k 


2 


2 

1 + sin$sin4>Q  + cos^cos^qCosAX 


(24) 


Calculation  of  the  Radius  of  the  Conformal  Sphere 

The  total  scale  factor  in  mapping  a point  from  the  ellipsoid  to 
the  common  coordinate  plane  is  the  product  of  k^  and  k^.  The  total 
scale  factor  is  expressed  as  follows. 


2Eccos<J> 


NcosL(l  + sin<j>sin$Q  + cos^cos^cosAX) 


(25) 


The  scale  factor  expresses  the  ratio  of  a linear  element  on  the 
common  coordinate  plane  to  the  corresponding  element  on  the  ellipsoid. 
The  scale  factor  therefore  represents  the  ratio  of  the  velocity  on 
the  common  plane  to  the  corresponding  velocity  on  the  ellipsoid.  A 
unity  scale  factor  is  highly  desirable  since  velocity  on  the  common 
plane  will  represent  actual  ground  speed.  The  scale  error  is  defined 
as  the  difference  between  the  scale  factor  at  a point  in  the  region 
and  a unity  scale  factor.  The  scale  error  e is  expressed  as  follows. 

2E  coa<P 

p a III  m .III  ■—  ■ — ■— ^ ^ ^ 

NcosL(l  + siniJ>sin$Q  + cos$cos$qCosAX) 

- E A - 1 (26) 

c 

The  extent  of  an  air  surveillance  region  is  defined  by  the  location 
of  the  radars  that  are  tied  into  it.  Since  the  scale  error  varies 
as  a function  of  the  location  and  separation  of  a point  from  the 
region  center,  it  is  highly  desirable  to  minimize  the  maximum  scale 


28 


errors  that  will  be  encountered.  From  equation  (26)  the  value  of 

the  scale  error  at  a point  can  be  varied  by  varying  the  radius  of 

the  conformal  sphere  Ec»  It  is  therefore  possible  to  obtain  both 

positive  and  negative  scale  errors.  To  minimize  the  magnitude  of 

the  largest  scale  error,  E is  chosen  such  that  the  magnitude  of  the 

c 

maximum  negative  scale  error  is  equal  to  the  maximum  positive  scale 
error.  The  radius  of  the  conformal  sphere  is  obtained  as  follows. 


"min 


+ £ 


max 


0 


EA  . -1  + EA  -1-0 

c min  c max 


E (A  . + A ) - 2 

c min  max 


E “t* 
c A 


. . + A 

min  max 


(27) 


where  A . and  A are  the  smallest  and  largest  A calculated  for 
min  max 

the  region  center  and  all  radars  tied  to  the  region.  Substituting 
equations  (19)  and  (24)  into  equation  (2*^)  the  value  of  A may  be 
calculated  as  follows. 


A 


(28) 


The  values  of  A and  A . for  a given  region  can  therefore  be 
max  min 

determined  from  Tables  II  and  III,  The  following  general  conclusions 
can  be  drawn  from  examination  of  the  tables. 

1,  Amin  usually  corresoonds  to  the  region  center  since  its 
value  of  kj  is  unity. 

2.  A usually  corresponds  to  the  most  distant  radar  since  it 

1RAX 

has  the  largest  angular  separation.  If  two  radars  have  the 

same  angular  separation,  A will  correspond  to  the  more 

max 

northerly  since  it  will  have  the  larger  kj^  value. 
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Table  IV  shows  values  of  E calculated  for  a region  center  at  a 

c 

geographic  latitude  of  45°  and  radar  sites  directly  north  and  south 
of  the  region  center.  Table  IV  shows  that  the  radius  of  the  conformal 
sphere  decreases  markedly  as  the  region  size  increases. 


Table  IV 

E Versus  Region  Size 
c 


Geographic  Latitude 
of  Most  Distant  Radar 

Distance  From 
Region  Center 

E 

c 

47° 

120.012 

3437.561 

49° 

240  ,10 

3435.787 

51° 

359.997 

3435.963 

53° 

479.970 

3429.088 

57° 

719.881 

3418.170 

61° 

959.744 

3403.004 

65* 

1199.563 

3383.553 

43° 

120.026 

3437.964 

41* 

240.064 

3436.596 

39° 

360.116 

3434.179 

37° 

480.179 

3430.714 

33" 

720.342 

3420.628 

29" 

960.553 

3406.314 

25" 

1200.808 

3387.749 

Effect  of  Ec/E 

s 

The  JSS  stereographic  ground  range  equation,  equation  (8),  and 
the  current  SAGE/BUIC  stereographic  ground  range  equation,  equation 
(9),  differ  by  the  scale  factor  Ec/Eg.  The  effect  of  E^E^  can  be 
expressed  as  the  difference  between  the  two  equations  as  follows. 
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Since  stereographic  ground  range  is  transformed  into  coordinates  on 
the  common  plane,  a difference  in  the  ground  ranges  will  result  in  a 
corresponding  misregistration  on  the  common  plane.  The  quantity  R 
in  equation  (29)  represents  an  approximation.  Since  the  design 
registration  error  budget  in  JSS  is  .18  nmi,  R must  differ  from  its 
actual  value  by  no  more  than  .18  nmi.  (A  detailed  review  of  the 
errors  induced  by  an  error  in  R is  beyond  the  scope  of  this  report.) 
Therefore  a difference  calculated  by  equation  (29)  greater  than  .36 
nmi  will  guarantee  an  unacceptable  registration  error  if  the  SAGE/ 
BUIC  stereographic  ground  r 2e  is  used.  To  show  the  effect  of 

E 

c/E  as  a function  of  region  size,  values  of  5 are  shown  in  Table 
s 

V for  the  radar  locations,  region  center  and  values  of  E£  indicated 
in  Table  IV.  The  value  of  R was  arbitrarily  chosen  to  be  100  and 
200  nmi.  Region  size  is  defined  as  the  distance  of  the  most  distant 
radar  from  the  region  center. 


Table  V 

Difference  Between  Ground  Range  Equations  8 and  9 versus  Region  Size 


Geographic  Latitude 
of  Most  Distant  Radar 

Approximate 
Region  Size 

E /F. 
c s 

6 (nmi) 
R**100nmi 

<S  (nmi) 
R“200nml 

47® 

120 

.9999 

.009 

.018 

49° 

240 

.9995 

.049 

.098 

51° 

360 

.9988 

.120 

.239 

53® 

480 

.9978 

.221 

.442 

57® 

720 

.9948 

.517 

1.033 

61® 

960 

.9906 

.937 

1.875 

65® 

1200 

.9852 

1,485 

2.969 

43® 

120 

.9998 

.021 

.042 

41® 

240 

.9993 

.072 

.145 

39® 

360 

.9985 

.154 

.308 

37® 

480 

.9973 

.266 

.533 

33® 

720 

.9942 

.581 

1.163 

29® 

960 

.9898 

1.018 

2.036 

25® 

1200 

.9842 

1.576 

3.152 
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Examination  of  Table  V reveals  that  for  R equal  to  200  nai, 
values  of  6 will  exceed  .36  nai  for  regions  sone where  between  360  nal 
and  480  nal.  This  Indicates  that  the  SACE/BUIC  stereographic  ground 
range  equation  will  produce  unacceptable  results  In  large  regions. 
Soae  of  the  current  SAGE/BUIC  regions  exceed  these  limits , and  all 
JSS  regions  will  exceed  these  limits  by  a considerable  margin. 
Section  III  of  this  report  is  devoted  to  depicting  the  registration 
errors  produced  by  the  SAGE/BUIC  and  JSS  stereographic  ground  range 
equations  In  the  seven  JSS  regions. 
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SECTION  III 


EFFECT  OF  THE  SAGE/BUIC  AND  JSS  STEREOORAPHIC  GROUND  RANGE 
EQUATIONS  ON  REGISTRATION 


INTRODUCTION 

Accurate  stereographic  projection  of  radar  data  onto  the  comon 
coordinate  plane  is  vital  to  the  operation  of  air  defense  and  air 
traffic  control  systems.  Large  registration  errors  seriously  down- 
grade the  performance  and  stability  of  the  active  tracking  algorithm. 
For  each  of  the  seven  JSS  regions,  simulated  radar  data  was  used  to 
demonstrate  the  effect  of  the  SAGE/BUIC  and  JSS  stereographic  ground 
range  equations  on  registration. 

SIMULATED  RADAR  DATA 

Radar  slant  range  and  azimuth  data  were  generated  for  four  air- 
craft locations  at  altitudes  of  30,000,  45, 000  and  60,000  ft  in  each 
JSS  region.  Appendix  1 indicates  the  algorithm  used  to  produce  the 
data.  For  a particular  radar,  a slant  range-azimuth  pair  was  gener- 
ated only  if  the  aircraft  was  within  250  nai  of  the  radar  and  above 
the  radar  horizon.  Aircraft  locations  were  chosen  so  that  slant 
ranges  for  the  reporting  radars  would  be  greater  than  170  nmi. 

Radar  slant  ranges  were  converted  into  stereographic  ground 
ranges  using  the  SAGE/BUIC  and  JSS  stereographic  ground  rsnge  e- 
quations  (equations  (8)  and  (9)) . Stereographic  ground  ranges  were 
transformed  into  coordinates  on  the  common  plane  using  equations  (16) 
and  (17).  Since  the  location  of  the  aircraft  was  known,  actual 
coordinates  on  the  common  plane  were  computed  using  equations  (13) 
and  (14).  The  registration  error  was  obtained  by  taking  the  magnitude 
of  the  difference  between  the  actual  coordinates  on  the  cowon  plane 
and  chose  obtained  by  the  SAGE/BUIC  and  JSS  stereograpblc  ground  range 
equations. 
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Appendix  II  presents  the  data  for  the  seven  JSR  regions . For 
each  region,  there  is  a Site  Data,  a Simulated  Radar  Data,  and  a 
Registration  Error  table.  In  addition  there  is  a figure  showing 
the  location  of  the  region  center,  radar  sites,  and  aircraft.  The 
following  information  is  given  in  the  tables. 

Site  Data  Table: 

1.  Approximate  geographic  latitude  and  longitude  of  the 
region  center  and  radar  sites. 

2.  Radius  of  the  earth  to  each  radar  site. 

3.  Coordinates  on  the  common  plane  for  each  radar  site. 

4.  The  radius  of  the  conformal  sphere. 

Simulated  Radar  Data  Table: 

1.  The  geographic  latitude  and  longitude  and  altitude  of  the 
aircraft  in  the  region. 

2.  The  designation  of  all  reporting  radars. 

3.  Slant  range  and  azimuth  data  for  all  reporting  radars. 

Registration  Error  Table: 

1,  Stereographic  ground  range  calculated  using  the  SAGE/BUIC 
and  JSS  ground  range  equations. 

2.  Coordinates  on  the  common  plane  obtained  from  the  SAGE/ 

BUtC  and  JSS  ground  range  equations. 

3*  Coordinates  on  the  common  plane  obtained  from  the  actual 
aircraft  locations. 

4.  The  registration  error  Induced  by  the  SAGE/BUIC  and  JSS 
ground  range  equations. 

Examination  of  the  tables  in  Appendix  IX  reveals  that  use  of  the 
SAGE/BUIC  stereographic  ground  range  equation  resulted  in  registration 
errors  that  exceeded  the  .18  tuai  JSS  registration  error  budget  in  all 
cases  tested.  Use  of  the  JSS  stereographic  ground  range  equation 
resulted  In  acceptable  registration  errors.  Table  VI  summarizes  the 
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worst  case  results  for  slant  ranges  of  approximately  180  nai  in  each 
of  the  seven  JSS  regions. 

The  omission  of  the  factor  E /E  is  responsible  for  the  large 

c s 

registration  errors  produced  by  the  SAGE/BUIC  stereographic  ground 
range  equation.  For  a particular  region,  the  value  of  Efi  is  constant. 
The  moat  southerly  radar  will  have  the  largest  value  of  Eg.  There- 
fore the  value  of  S as  calculated  by  equation  (29)  will  be  largest 
for  the  most  southerly  radar.  This  indicates  that  the  worst  case 
errors  will  be  produced  by  the  most  southerly  radar,  and  the  best 
case  errors  by  the  most  northerly.  Further  examination  of  Appendix 
XI  supports  this  conclusion. 
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SECTION  IV 


CONCLUSIONS 

The  scale  factor  E /E  should  be  included  in  the  stereographic 
c s 

ground  range  equation  to  avoid  large  registration  errors.  This  is 
especially  important  in  large  regions.  The  JSS  stereographic  ground 
range  equation  for  processing  returns  with  height  data  is: 


The  present  SAGE/BUIC  stereographic  ground  range  equations 
should  be  modified  to  reflect  the  scale  factor  Sc/Es.  The  modifi- 
cation would  require  a change  in  the  adaptation  parameters  on  a site 
by-site  basis. 
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Preceding  page  blank 


APPENDIX  I 


GENERATION  OF  SLANT  RANGE  AND  AZIMUTH  DATA 


Radar  slant  range  and  azimuth  data  are  calculated  for  aircraft 


locations  using  vector  operations.  The  vector  V from  the  center 
of  the  earth  to  a point  on  or  above  the  ellipsoid  is  calculated  as 


follows. 


where : 


X « 
v 

Z - 

V 


(N  + H)cos  LcosX 
(N  + H)cosLsinX 
{N(l  - e2)  + HjsinL 


(A-l) 


H is  the  height  of  the  point  above  the  ellipsoid 
X is  measured  positive  east  of  the  prime  meridian 
Hie  slant  range  S is  computed  as  follows. 

S - T - R 

S - j£|  (A-2) 

where; 

T,  R are  vectors  from  the  center  of  the  earth  to  the  aircraft 

tiw. 

and  radar  respectively. 

The  aircraft  must:  be  above  the  radar  horizon  for  a radar  return  to 
be  possible.  The  aircraft  is  above  the  radar  horizon  if; 

£ • Z » 0 

where: 

Z is  a unit  vector  directed  along  the  zenith  of  the  reporting 
radar  as  follows; 


Pncotaf  pan  Mask 
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rX  » cosL  cosA 
- r r 


Z ■•{Y  ■ cosL  sinA 
z r r 


Z * sinL 
z r 


(A-3) 


Lr»  *r  are  the  geographic  latitude  and  longitude  of  the  report- 
ing radar 

The  azimuth  angle  6 is  calculated  as  follows* 

S‘El 


0 • tan 


-1 


S»N 


(A-4) 


where: 


E is  a unit  vector  directed  due  east  of  the  reporting  radar  as 
follows. 


fX_  - -sinA 
E r 


Y_  ■ cosA 
E r 


(A-5) 


M is  a unit  vector  directed  due  north  of  the  reporting  radar  as 
follows . 


- -sinL^cosA^ 


N -< 


Y ■ -sinL  sinA 
H r r 


cosL 


(A-6) 
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APPENDIX  II 


SIMULATED  DATA  FOR  THE  SEVEN  JSS  REGIONS 

Radar  slant  range  and  azimuth  data  were  generated  for  four  air- 
craft locations  at  altitudes  of  30,000,  45,000  and  60,000  feet  in 
each  JSS  region.  For  a particular  radar,  a slant  range  azimuth  pair 
was  generated  if  the  aircraft  was  within  250  nmi  of  the  radar  and 
above  the  radar  horizon.  Radar  slant  ranges  were  converted  into 
stereographic  ground  ranges  using  the  SAGE/BUIC  and  JSS  ground  range 
equations.  Stereographic  ground  ranges  were  transformed  into  co- 
ordinates on  the  common  coordinate  plane.  The  actual  coordinates  on 
the  conmton  plane  were  also  obtained.  The  registration  errors  Induced 
by  the  SAGE/BUIC  and  JSS  ground  range  equations  were  calculated. 

For  each  JSS  region,  there  is  a Site  Data,  Simulated  Radar  Data, 
and  Registration  Errors  table.  Tn  addition,  there  is  a figure 
showing  radar  site,  aircraft,  and  region  center  locations. 


Explanation  of  Tables 

The  Site  Data  Table  shows  the  following. 

a)  Approximate  latitude  and  longitude  of  the  radar  sites  and 

region  center. 

b)  The  earth  radius  to  each  radar  site  calculated  using 

equation  (12). 

c)  Coordinates  on  the  common  plane  for  the  radar  sites  calcu- 

lated using  equations  (13)  and  (14). 

d)  The  radius  of  the  conformal  sphere  calculated  using  equation 

(27). 

The  Simulated  Radar  Data  Table  shows  the  following. 

a)  The  latitude,  longitude  and  altitude  of  the  aircraft. 

b)  The  slant  range  and  azimuth  for  all  reporting  radars.  A 
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radar  la  a reporting  radar  if  the  target  is  within  250 
nml  and  above  the  radar  horizon.  Slant  range  and  azimuth 
data  are  calculated  using  the  algorithm  In  Appendix  I. 

The  Kegistratlon  Error  Table  shows  the  following. 

a)  The  ntereographlc  ground  range  calculated  using  the  SAGE/ 

BOIC  and  JSS  ground  range  equations  (equations  (8)  and 

(9)). 

b)  Coordinates  on  the  common  plane  obtained  from  the  SAGE/BUIC 

and  JSS  stereographic  ground  range  equations  using 
equations  (16)  and  (17). 

c)  Coordinates  on  the  common  plane  obtained  from  the  actual 

location  of  the  aircraft  using  equations  (13)  and  (14). 

d)  The  registration  error  Induced  by  the  SAGE/BUIC  and  JSS 

stereographic  ground  range  equations. 
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